Mon. Not. R. Astron. Soc. 000, [HS] (2007) Printed 22 May 2008 (MN WI^ style file v2.2) 



Exploring the photometric signatures of magnetospheres 
around Helium-strong stars 



oo R. H. D. Townsend* 

Bartol Research Institute, Department of Physics & Astronomy, University of Delaware, Newark, DE 19716, USA 

o 



<N 

6 



(N 
> 

rn 

o 

oo 

o 



- 

X 



Received: Accepted: 



ABSTRACT 

The photometric variations due to magnetically confined material around He-strong 
stars are investigated vifithin the fi-amework of the Rigidly Rotating Magnetosphere 
(RRM) model. For dipole field topologies, the model is used to explore how the mor- 
phology of light curves evolves in response to changes to the observer inclination, 
magnetic obliquity, rotation rate and optical depth. 

The general result is that double-minimum light curves arise when the obliquity 
and/or inclination are close to 90°; no light variations are seen in the opposite limit; 
and for intermediate cases, single-minimum light curves occur. These findings are 
interpreted with the aid of a simple, analytical torus model, paving the way for the 
development of new photometric-based constraints on the fundamental parameters of 
He-strong stars. Illustrative applications to five stars in the class are presented. 

Key words: stars: chemically peculiar - stars: fundamental parameters - stars: 
magnetic fields - stars: rotation - stars: variables: other ~ techniques: photometric 



1 INTRODUCTION 

The Helium-strong stars are small class of chemically 
peculiar B stars, interpreted as an extension of the 
Ap/Bp phenomen on to higher effective temperatures 
(|Osmer fc Potorsg^ 1 1974 ). They are characterized by el- 
evated and spatially inhomogeneous photospheric Helium 
abundances, and in most cases al s o exh ibit strong mag- 
netic fields (see IBohlender et al.l 1 19871 . and references 
therein). Spectroscopy of their time- varying Helium ab- 
sorption lines often reveals that the surface abundance 
distribution is strongly corr e lated with field topology 
IBohlender fc Landstreetlll988l : IVetdll990l '). These correla- 
tions lend support to models for the abundance anomalies 
that rely on the interplay between radiative and magnetic 
forces, ei ther i n the framework of elemental diffusion (e.g., 
Michau d'ITgTOh. o r in t he context of wind mass loss (e.g.. 
Hunger fc GrootdligggD . 

Of the He-strong stars exh ibiting photometric variabil- 
ity l|Pedersen fc ThomserJ 0*9771 ). at least some can success- 
fully be explained by invoking rotational modulation of the 
visible a bundance inhomogen eities - so-called 'spot' models 
(see, e.g.. lKrticka et al]|2007l ). However, a different approach 
is required for other objects in the class, most notably the 
B2pe star a Ori E. The light curve of this star has a d ouble- 
eclipse morphology l|Hesser. Ugarte fc Moreiiol 19771), l ead- 
ing some authors (e.g., Hesser. Walborn fc Ugart^l 19761 ) ini- 
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tially to suggest that the star is a mass-transfer binary. In 
fact, as further observational data have become available, 
an empirical picture for the star's ph otometr ic and spec- 
troscopic variability has emerged (Groote fc Hunger. 198 j). 
that augments a spot model with a pair of circumstellar 
clouds. These clouds absorb continuum photons as they 
transit the stellar disk, giving rise to the distinctive photo- 
metric minima; however, once off the disk they reveal them- 
selves in Hq line profiles as blue- and red-shifted emissio n 
components (e.g.. IWalborn fc Hesser|[l976l : |PederserJ[l979l ). 

Recently. iTownsen d fc Owockil (|2005l . hereafter TO-05) 
have developed a Rigidly Rotating Magnetosphere (RRM) 
model that provides the theoretical basis for this empirical 
picture. The model conjectures that plasma in the star's ra- 
dial ively driven wind, channeled by nearly rigid field lines, 
tends to accumulate at points where the effective (gravita- 
tional plus centrifugal) potential is at a local minimum. This 
steady accumulation leads to the development of a magneto- 
sphere resembling a warped disk, which is supported by the 
centrifugal force and co-rotates rigidly with the star. For 
an oblique dipole field topology, the RRM model predicts 
that the magnetospheric density is highest in two regions 
- 'clouds' - situated along the twin intersections between 
the magnetic and rotational equators. In the specific case of 
a Ori E, the photometric and Ha variations synthesized us- 
ing the mo del show very good agreement with observations 
of the star (|Townsend. Owocki fc Grootell2005l ). 

In the present paper I apply the RRM model in a 
broader context, by exploring how the photometric signa- 
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tures of magnetospheric plasma depend on fundamental 
parameters such as the observer inclination and magnetic 
obliquity (^. The results from this parameter-space explo- 
ration are interpreted with the aid of a simple, analytic torus 
model (ill), and I discuss how they might be applied to con- 
strain the parameters of He-strong stars O. The principal 
findings are then summarized (^Sl)- 



2 EXPLORATION OF PARAMETER SPACE 
2.1 Light-curve synthesis 

The RRM model as applied to a dipole magnetic topol- 
ogy is used to synthesize the light curves used for ex- 
ploring parameter spa ce. The same procedure described by 
I Townsend et al.l (|2005l . their §3.2) is followed; note therefore 
that the variability arising from abundance inhomogeneities 
is not considered. It is assumed that there is no decentering 
of the magnetic field, a = (eqn. 1, ibid.), and a canonical 
energy ratio parameter e* — 10~'^ (cf. TO-05, eqn. 38) is 
adopted throughout. 

The remaining free parameters of the model are then 
the inclination i; the magnetic obliquity /?; the ratio 



(1) 



between the rotation frequency SI and the critical frequency 
f2c; and an optical depth scale 



To = 



K M t U) 



4/3 
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This latter expression groups the opacity k (assumed con- 
stant) with various terms appearing in the RRM prescrip- 
tion for the density (TO-05, eqn. 35). The additional oj*''"^ 
term comes from considering the rotational dependence of 
the total magnetospheric mass column between star and ob- 
server; its inclusion ensures that models having the same 
tq but difi^ering uj exhibit approximately the same degree of 
light variations. 



2.2 Varying /3 and i 

To begin the exploration of parameter space, a grid of light 
curves is synthesized in 5° increments over the full intervals 
0° ^ i ^ 90° and 0° ^ /3 90°. In each case an interme- 
diate rotation rate oi uj = 0.5 is assumed, and an optical 
depth scale ro = 10 is chosen to produce variability on the 
~ 10% level typical to He-strong stars. Figure [T] illustrates 
the resulting light curves, plotting the flux F (in units of the 
unobscured stellar flux Ft) against the rotational phase <j!0. 
For clarity only the curves at even multiples of 5° are shown, 
and the i < 30° cases are omitted because none exhibit any 
signiflcant variability. 

The flgure clearly demonstrates how the character of the 
light curves is altered as i and f3 are changed. In the limit 
where either of these parameters is large, the curves exhibit a 



^ In the present work, </> = 0° corresponds to magnetic maxi- 
mum; this differs from the convention adopted by 'Townsend et al.| 
1I2OO5I ). who followed historical precedent (e.g., Hesscr et al., 1977i) 
by choosing the primary light minimum as phase zero. 



clear double-minimum (hereafter, '2-m') morphology. In the 
opposite limit, and in particular for i < 30° as already men- 
tioned, a non-varying, constant morphology ('0-m') is seen. 
For intermediate cases, falling approximately along a band 
extending from {i,/3} = {30°, 90°} to {i,/3} = {90°, 0°}, a 
single light minimum ('1-m') occurs at phase <j!> = 180°. 

To understand the origin of these three morphologies. 
Fig. [2] shows maps of the column density D in the near-star 
regions of the circumstellar environment, for three conflgu- 
rations having /3 — 50° and i = 30° , 50° , and 90° . The maps 
are calculated from an expression of the form 



D{xo,yo) 



p(Xo, yo,Zo) dZo, 



(3) 



where (a;o, yo, Zo) are Cartesian coordinates in the observer's 
reference frame, with Zo directed toward the observer, and 
p is the local density predicted by the RRM model (cf. TO- 
05, eqn. 35). Panel (i) of the figure, plotting the i = 30° 
case, is typical of low-inclination systems. With the observer 
situated close to the rotation pole, the plasma trapped in 
the magnetosphere - lying in a disk between the rotational 
and magnetic equatorial planes - has no opportunity to pass 
between the observer and the star. Hence, an unvarying 0-m 
light curve is seen in these systems. 

For panel (iii) the converse is true: with the observer in 
the rotational equatorial plane, the magnetospheric plasma 
transits the star twice per rotation cycle, and thus a 2-m 
light curve arises. In this particular case, symmetry requires 
that the minima occur at phases (p — 90° and 270°, with a 
half-cycle difference between them. However, as discussed in 
further detail below, a smaller phase difference A(j) is found 
for 2-m configurations having i < 90° and /? < 90°. 

Panel (ii) illustrates the intermediate, i — 50° case, 
where the observer is not sufficiently close to the rotational 
equator for magnetospheric plasma to pass fully in front of 
the star. Instead, when the magnetic field is pointed away 
from the observer the plasma grazes the lower portions of 
the stellar limb, causing a 1-m light curve with the minimum 
centered at phase (f) = 180°. This light curve is character- 
ized by a reduced overall absorption, relative to the 2-m case 
shown in panel (iii). 

As can be seen from the l3 = 50° row of light curves in 
Fig.[T] a smooth evolutionary sequence links the cases shown 
in panels (i), (ii) and (iii) of Fig. [5] As the observer moves 
from the rotational equator toward smaller i, the depth 
(= F, — F) at the twin light minima of the (initially, 2-m) 
light curve decreases; whereas, the depth at the interven- 
ing light maximum at (f) = 0.5 increases. Eventually, these 
depth values cross over, with the result that the light curve 
transitions to a 1-m morphology. Following this switch-over 
(which in the /3 = 50° case occurs at i ~ 60°), the depth at 
the single remaining minimum at = 180° declines toward 
lower i, and eventually vanishes in a further transition to a 

0- m morphology. 

The progressive distortion of light curves, as the 2-m to 

1- m transition is approached, has an infiuence on the phase 
difference between the minima. To explore this infiuence, 
each light curve in the grid is analyzed using a simple 
algorithm that determines its morphology, and - if 2-m - 
measures the phase difference. The algorithm involves three 
steps: 
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Figure 1. Light curves synthesized with varying inclination i and obliquity /3, in each case plotting the flux F (in units of the unobscured 
stellar flux F*) as a function of the rotation phase 0; a rotation rate u) = 0.5 and an optical depth scale tq = 10 are assumed. The dotted 
horizontal lines indicate the level F = F». 
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(i) I = 30° (ii) i = 50° (iii) i = 90° 




-4 -2 2 4 -4 -2 2 4 -4 -2 2 4 

Figure 2. Maps of the column density D in the observer's image plane (xo,yo), for three RRM configurations: (i) {i, /?} = {30° , 50°} at 
a rotational phase = 81°, (ii) {i,f3} = {50°, 50°} at </> = 171°, and (iii) {i,/3} = {90°, 50°} at (f> = 81°. (The phases are ofl'sct slightly 
from the 'round-number' values 90° and 180°, to improve visual clues to the 3-D geometry). In all three cases, lo = 0.5. Black/opaque 
corresponds to highest column density, and white/transparent to lowest, with intermediate levels shown in orange. The central star is 
coloured blue, and the axes are drawn as arrows — yellow/upright for the rotation axis, and magenta/oblique for the magnetic axis. 
Moreover, selected field lines having a summit radius of 5 Rp (where Rp is the stellar radius at the rotational poles) are indicated in 
green. 

A0 



45° 90° 135° 180° 




Figure 3. The phase difference A<j) between the minima in 2-m light curves, plotted in the plane for three differing rotation rates; 
in every case, tq = 10. The curve in each panel shows the contour A0 = 150° , discussed in i|4.2.2l Regions of the plane associated with 
0-m and 1-m curves are shown in monochrome and labeled as '0' and '1', respectively. 

(i) If the total flux variation AF = max(_F) — min(_F') is both satisfy F — min(F) > 0.15A_F, then the Ught curve is 
less than 0.01 , then the Ught curve is classified as 0-m. classified as 2-m. 

(iii) Otherwise, the hght curve is classified as 1-m. 



(ii) Otherwise, if the hght curve exhibits two fiux min- A minor complication arises when minima are cusped 

ima, and these minima are separated by flux maxima that (see ^2.3l and Fig.|4]); then, the adopted phase of a minimum 
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Figure 4. Light curves for a {«,/?} = {70°, 50°} configuration, 
synthesized with tq = 10 and rotation rates ui = 0.25, 0.50, and 
0.75. For clarity, the w = 0.50 curve is displaced downward by 
0.02F*, and the uj = 0.75 curve by 0.04F,. 



Figure 5. Light curves for a = {70°, 50°} configuration, 

synthesized for u> = 0.50 and optical depth scales tq = 5, 10, 
and 15. For clarity, the tq = 10 curve is displaced downward by 
0.02 F,, and the tq = 15 curve by 0.04 F*. 



is determined by averaging the phases of the local minima 
on either side of the cusp. 

Figure |3] plots the results of this analysis, showing the 
Aff) values for the 2-m light curves in the plane, and also 
indicating where in the plane 1-m and 0-m curves occur. 
(The middle panel of the figure shows the ui = 0.5 case con- 
sidered in this section; the other two panels, corresponding 
to a; = 0.25 and w = 0.75, are discussed in the following 
section). As before, the i < 30° region of the plane is uni- 
formly 0-m and is therefore not shown; but in the present 
case, the /3 < 20° region is also omitted, because there the 
analysis algorithm gives unreliable results. This is more due 
to the difficulty in sensibly classifying the light curves in 
this region (in particular, see the {i,/3} ~ (80°, 10°} and 
{90° , 10° } cases shown in Fig. [TJ , than to any particular 
limitation of the algorithm. 

The figure confirms the previous finding that the 0-m 
and 2-m light curves are situated toward small and large 
angles, respectively, with the 1-m curves falling in between. 
However, it also reveals that A(p depends quite regularly on 
i and (3. When either of these parameters is equal to 90°, 
the phase difi'erence is exactly one-half of a cycle; as already 
mentioned, this is required by symmetry. Away from these 
special cases, and toward the boundary between the 2-m and 
1-m regions, A(j) is constant along contours running from 
low-i/high-/3 corner down to the high-i/low-/3 corner. The 
origin of this behaviour is explored further in ^ 



2.3 Varying cu 

The three panels of Fig.O showing the light-curve morpholo- 
gies in the j-/3 plane for ui — 0.25, 0.5, and 0.75, illustrate 
some of the effects of varying the rotation rate. The most 
obvious trend is the expansion of the 1-m region as uj is in- 
creased. This comes almost entirely at the expense of the 
0-m region; almost no change is seen in the 2-m region, ei- 
ther in terms of its extent, or in terms of the Acp values 



shown. Again, the reasons for such behaviour are explored 
further in iJJ] 

Focusing now on an individual case. Fig. |4] plots the 
light curves for a {i,l3} — {70°, 50°} configuration [chosen 
as half-way between the cases shown in panels (ii) and (iii) 
of Fig. [2], synthesized for the same three uj values above. As 
already noted for Fig. [S] the rotation rate has only a minor 
impact on the phase difi'erence between the light minima; 
the measured differences are — 118° for uj — 0.25, A0 — 
116° for UJ = 0.50, and Aij) = 124° for uj = 0.75, amounting 
to only a 10% variation from a factor-three increase in uj. 

However, a significant increase occurs in the width of the 
individual light minima toward more-rapid rotation. This in- 
dicates that the magnetospheric plasma is taking relatively 
longer to transit across the star, and follows from the ap- 
proximate scaling of the RRM inner-edge radius with the 
Kepler radius. 



(cf. TO-05, eqn. 12). As uj increases, both Rk and the inner 
edge moves in toward the star. Thus, the plasma takes a 
larger fraction of a rotation cycle to transit, and broader 
light minima ensue. 

An additional effect of rapid rotation, seen in the uj = 
0.75 panel of Fig. |4l is to produce cusping at the minima of 
the light curve. This effect, which is the characteristic hall- 
mark of an optically thick yet geometrically thin disk, arises 
because the scale height hm of the disk-like magnetosphere 
varies as ~ R^l'^ ~ UJ-'' (cf. TO-05, eqn. 29). 

2.4 Varying tq 

For the sake of completeness Fig. [5] plots the light curves 
for the same {i,/3} — {70°, 50°} configuration considered 
previously, but synthesized this time with a fixed rotation 
rate uj — 0.50 and varying optical depth scale: ro = 5, 10 
and 15. The maximum depths of the light curves - 0.037 F,, 
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Figure 6. The torus model for a {i,vt,Rt} = {70° , 50° , 3.5 Rp} 
configuration, seen at rotational phase = 92°; the torus is col- 
ored orange, while the star and axes are as in Fig. [2] The dotted 
line shows the ellipse corresponding to the sky projection of the 
torus; this ellipse is centered at point O, and intersects the stellar 
limb at points A and B. The arc AB has length £, and the line 
OM - where M is the midpoint of the arc - is the semi-minor axis 
b of the ellipse. 



0.069 and 0.099 , respectively - increase approximately 
linearly with to- This result is a natural consequence of the 
fact that the magnetospheric densities (cf. TO-05, eqn. 35) 
are proportional to tq. As in the preceding section, however, 
the change in phase difference is small: in the direction of 
increasing to, Aij) = 120°, 116° and 114°. 



3 A TORUS MODEL 

The previous sections explore how RRM light-curve mor- 
phology is sensitive to parameters such as i, j3, and uj. A 
number of interesting dependencies emerge, as seen for in- 
stance in Fig. [3l To gain a measure of insight into these 
dependencies, the present section introduces an idealized, 
analytic torus model for the formation of the light curves. 

In this model, flux variations arise due to obscuration 
of a spherical star by an optically thick, oblique, co-rotating 
torus of circumstellar material. Assuming that the torus is 
narrow compared to the stellar radius, the observed flux at 
any instant can be approximated by considering the fraction 
of the stellar disk obscured: 



(5) 



Here, i?* is the stellar radius, dt R*) is the torus cross- 
sectional diameter, and I is the length of the torus arc (if 
any) covering the stellar disk. For Rt not too close to the 
stellar radius _R*, this arc length is approximated by 



2y'max(i?2 _fo2 0), 



(6) 



where h is the semi-major axis of the ellipse that is the pro- 
jection of the torus on the plane of the sky (see Fig. [6]). If 



the torus plane is tilted by an angle Vt to the rotation axis, 
it can be shown that 



b = Rt \ cos i cos Vt + sin i sin vt cos ( 



(7) 



where Rt is the radius of the torus. 

Even before light curves are synthesized, these expres- 
sions can be used to demonstrate that the torus model re- 
produces the same three light-curve morphologies found pre- 
viously. For any light variations whatsoever to occur, 6 must 
be smaller than i?, at some point during the rotation cycle. 
Assuming that 0° ^ i ^ 90° and Q° ^ Vt ^ 90°, this is is 
equivalent to the requirement 



cos(i + vt)< 



R^ 
Rt' 



(8) 



If this condition is fulfilled, then the light curve exhibits 
a single or double minimum morphology depending on 
whether the torus passes across the centre of the stellar disk: 



cos(i -I- Vt) 



> 
< 



1- m 

2- m 



(9) 



In the 1-m case, the fiux minimum always occurs at phase 
(h — 180°. In the 2-m case, the minima arise when 



cos i cos Vt — ~ sin i sin Vt cos i 
so that the phase difference is given by 



(10) 



A(fi — 2cos (coticotut) (11) 

A sim ilar expression has been derived by IShore fc Brownl 
(|l990l . their eqn. 2), but in their case f3 appears in the place 
of Vt. 

To now relate the foregoing analysis back to the RRM 
model, representative values are chosen for the torus radius 
Rt and tilt angle Vt. For the former, the Kepler radius Rk 
(cf. eqn. is the appropriate choice; for the latter, Ap- 
pendix (IX)) argues that 



V = (5 — tan 



sin 2/3 
5 + cos 2/3 



(12) 



furnishes a good characterization of the mean tilt of the 
magnetospheric disk. Adopting these values, and further as- 
suming a diameter dt = O.li?, , eqns. ([5HZ| are used to syn- 
thesized a. 13 — i — UJ grid of light curves. These data are 
analyzed using the algorithm described in i]2.2l and Fig. [7] 
plots the results in the same format as Fig. [3] 

The two figures shown a good measure of qualitative 
agreement. Specific characteristics of the RRM model that 
are correctly reproduced by the torus model include: 

• the three differing light-curve morphologies (0-m, 1-m 
and 2-m), and their general locations in the {i, /?} plane; 

• the fact that the phase differences of 2-m light curves 
are insensitive to rotation rate; this follows from eqn. (|lip . 
which has no dependence (explicit or implicit) on uj\ 

• likewise, the fact that the boundary between the 1-m 
and 2-m regions is insensitive to the rotation rate; this fol- 
lows from eqn. (O, which again has no dependence on u}\ 

• conversely, the fact that the boundary between the 0-m 
and 1-m regions does depend on the rotation rate, moving 
toward lower i as w is increased; this is due to the identifi- 
cation Rt = i?K, which introduces an implicit dependence 
on Lo into eqn. ([8]); 
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CJ = 0.5 



CJ = 0.75 




Figure 7. As for Fig. |3] except that the hght curves are calculated using the torus model rather than the RRM model; in every case, 
dt = 0.1 R,. 



• finally, although not shown in the figures, the increase in 
the width of 2-m light minima toward more-rapid rotation. 

At a quantitative level, the agreement between the two 
models is more limited; for instance, at large /3 the torus 
model predicts an extension of the 2-m region to low inclina- 
tions, which is not seen in the RRM model. [In this particular 
case, the discrepancy arises because high-/? magnetospheres 
are punctured by cone-shape voids over the rotational poles 
(see TO-05, fig. 3), and therefore do not produce light varia- 
tions when viewed from low inclinations.] Such discrepancies 
mean that the torus model should not be used as the sole 
basis for data analysis. However, the buUeted list above is 
ample proof that - as a tool for data interpretation - the 
model has many useful insights to offer. 



4 DISCUSSION 

4.1 Photometric diagnostics 

The preceding sections demonstrate that the light variations 
arising from a RRM depend on fundamental stellar param- 
eters in a way that is simple, regular and straightforward to 
understand. Given these qualities, it is reasonable to conjec- 
ture that observations of these variations, in He-strong stars 
and similar objects, could be used in tandem with other di- 
agnostics such as magnetic field measurements to establish 
constraints on the parameters of these stars. 

The success of this endeavour rests to a large part 
on the validity of the RRM model, since it is from dia- 
grams such as Fig. [3] - constructed with the aid of the 
model - that reliable diagnostics are developed. In this 
respect, the close correspondence between the observed 



and s ynthesized Ha variations of a Ori E (|Townsend et al.l 
200^ has lent considerable empirical support to the model. 
Moreover, recent rigid-field hy drodynamical simulations by 
iTownsend. Owocki fc ud-Doul a (2007 ), and magnetohydro- 
dynamical simulations by ud-Doula. Owocki fc Townsendl 
(2007), have provided confirmation of the theoretical un- 
derpinnings of the model. 

Of equal importance, however, is the problem of disen- 
tangling the light variations arising from photospheric abun- 
dance inhomogeneities on the one hand, and from magne- 
tospheric plasma on the other. For certain configurations, 
these two mechanisms can produce variability that appears 
ostensibly the same. For instance, the {i,/3} = {70°, 20°} 
magnetospheric light curve in Fig.[T] showing a single, broad 
minimum, might easily be confused with the rotational mod- 
ulation arising from a spot-like photospheric inhomogene- 
ity. In situations like these, other diagnostics such as spec- 
troscopy or polarimetry must be used to lift the degeneracy 
and correctly identify which mechanism is responsible for 
the observed variability. 

In the majority of cases, however, the morphologies of 
the light curves produced by the alternative mechanisms 
are rather different. Especially toward large i and/or f3, 
the light variations caused by magnetospheric plasma tend 
to be much more rapid (relative to the rotation period) 
than could plausibly be achieved by photospheric inhomo- 
geneities. Thus, as a simple rule, faster variations can be at- 
tributed to the magnetosphere, and slower variations to the 
photosphere. By way of illustration, the following section 
applies this heuristic to five He-strong stars, in some cases 
examining whether the constraints provided by Fig. [3] are 
consistent with literature values for these stars' parameters. 
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and in other cases using such constraints to make predic- 
tions. 

4.2 Application to He-strong stars 

4.2.1 HD 37017 

The hght variations of this star (e.g., iBolton et al.l [iQOSi ) 
show a smooth, sinusoidal character, fav ouring a photo- 
spheri c rather than magnetospheric origin. iBohlender et ahl 
(|l983) derive 42° < /3 < 59° and 23° < i < 37° from mag- 
netic and period measurements; with reference to Fig. [S] 
these values mostly fall in the 0-m region and are there- 
fore consistent with the apparent lack of a magnetospheric 
signature. 

4.2.2 HD 37479 (a Ori E) 

As already discussed in ^ this star exhibits distinctive 2-m 
light variations that are consistent with a magnetospheric 
origin. That said, certain aspects of its light curve, such as 
the unequal depths of the minima and the 'emission' feature 
after the secondary minimum, appear more likely be of pho- 
tospheric origin (see Towns end et al.ll2005l : note that these 
authors explained the depth asymmetry by instead invoking 
a decentered dipole). 

Nevertheless, this photospheric signal is unlikely to have 
much effect on the measured phase difference /S.<j) ~ 150° 
between the two minima. This value is indicated in each 
panel of Fig. [3] by the contour curves. Along these con- 
tours the combine d angle 15 -\- i var i es be tween ~ 105° and 
~ 135°, which led iTownsend et al.l (|2005l ) to adopt the ap- 
proximate relation /3 -I- i « 130° in developing their model 
for the star. This constraint is cons istent with the p aram - 
eters « {70°, 56°} derived bv lShore fc Brownl (I199G| ) 

from consideration of magnetic and lUE data. 

4.2.3 HD 37776 

This star has a rather non-sinusoidal light curve, with a sin- 
gle, very broad minimu m suggestive of a phot ospheric origin. 
Detailed modeling by iKrticka et al.l l|2007l ) demonstrates 
that, indeed, the photosphe ric abundance dis t ribut ion deter- 
mined spectroscopically bv lKhokhlova et al.1 (|200G| ) can also 
faithfully reproduce the star's light curve. Note that the star 
has a quadrapolar field topology (|Thompson fc Landstree^ 
1 19851 ): although the RRM model can in principle be applied 
to any arbitrary field topology, studies so far (including the 
present paper) focus on the simplest and commonest case of 
a dipole field. Thus, the absence of a clear magnetospheric 
signature in the light curve of HD 37776 cannot, without 
further modeling, be used to constrain i or /3. 

4.2.4 HD 64740 

This star is the brightest in the He-strong class, 
yet it exhibits almost no photometric v ari ations 
llPedersen fc ThomsenI 1 19771 ). IBohlender erall (Il987l ) 
derive /? < 76° and i > 41° from magnetic and period 
measurements. Given the absence of any variability, the 
boundary between the O-m and 1-m regions in Fig. [3] 



(assuming uj = 0.5) can be used to further constrain the 
parameters as /3 + 2i < 150° . 

4.2.5 HD 182180 

This star has recently been shown to be He- strong, and to 
exhibit Ha and photometric variations, by iRivinius et al.l 
iiool). Observational data remain scarce, but the 2-m, 
A<^ = 180° morphology of the double-wave light curve is 
consistent with P — 90° and/or i = 90°. The light min- 
ima are rather too broad to convincingly rule out a photo- 
spheric origin, but given the star's extreme rotation period 
of 0.521 d, the widths of the minima are wholly consistent 
with a near-critical, ~ 1 configuration. 



5 SUMMARY 

The RRM model has been used to explore how the light 
variations due to magnetospheric plasma depend on the pa- 
rameters /3, i, UJ, and tq (!J21)- Double-minimum morpholo- 
gies are found toward large values of the inclination i and/or 
obliquity (3; single-minimum and non-varying morphologies 
occur for progressively smaller values of these parameters. In 
the double- minimum cases, the phase difference Acp between 
the minima depends primarily on i and /3, with almost no 
sensitivity to the rotation rate lu or optical depth scale tq. 

Insight into these findings has been developed with the 
aid of an analytic torus model (Q. Although not in com- 
plete quantitative agreement with the RRM model, the torus 
model is able to explain all of the qualitative results from the 
parameter-space exploration. Applying these results, con- 
straints on the obliquity and inclination of four He-strong 
stars are derived from their light curves (311); these are found 
to be generally consistent with values quoted in the litera- 
ture. 
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For field lines extending far outside the Kepler radius (i.e.. 
7^1), the solution is then found as 



sin 2/3 
5 + cos 2/3 ' 



(A4) 



The local tilt angle v of the disk with respect to the magnetic 
axis is determined from 



tanv — ~ — 



sin 2(3 
5 + cos 2/3 ■ 



(A5) 



Thus, the tilt angle with respect to the rotational axis fol- 
lows as 



f3 + V = (3 — tan 



sin 2/3 
5 + cos 2(3 



(A6) 



Strictly speaking, this value cannot be applied across the 
whole of an RRM disk, due to the warping that arises when 
/3 > 0°. However, as can be seen from Fig. 3 of TO-05, the 
disk tilt in the = 90° , 270° meridional planes furnishes a 
good approximation to the mean disk tilt, and therefore is 
well suited to characterizing the overall disk orientation. 



APPENDIX A: THE DISK TILT ANGLE 

As discussed in SjlJ the RRM model predicts magnetospheres 
that resemble a warped disk. A characteristic disk tilt angle 
V with respect to the rotation axis can be obtained by con- 
sidering the locus formed by the minima of the dimensionless 
effective potential. 



1 7^ sin'' ( 



7 sm 



(sin (3 



sin^ 6 sin^ th + 



cos 9 + cos (3 sin 9 cos ( 



(Al) 



(cf. TO-05, eqn. 22). Here, 8 is the magnetic colatitude on 
the dipole field line having summit radius 7 (in units of the 
Kepler radius Rk; cf. eqn. 14)1 and magnetic azimuth The 
conditions for a potential minimum are that 



0, 



> 0. 



(A2) 



As discussed by TO-05, algebraic solution of these equations 
for general (3, 7, and (p is not possible. However, the (p = 
90°, 270° meridional planes are exceptions - for these, it is 
straightforward to demonstrate that minima occur at the 
magnetic equator, 9 = 90°. 

Immediately adjacent to these special cases, for instance 
in the (p — 90° -I- 5^ meridional plane for \5^\ <C 1, the 
continuity of the disk requires that the potential minimum 
should be at S = 90° + So, where ^ 1 also. Substituting 
this trial solution into the zero-gradient condition in (|A2|) . 
and neglecting terms of quadratic- or higher-order in 5^ and 



leads to the equation 



7^ sin 2(3 



5s + 



2 7^ (5 + cos 2/3) 
"7+ 2 



0. (A3) 



^ Not to be confused with the rotational phase 4>. 



